Relaxation of frustration and gap enhancement by the lattice distortion in the A chain 
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We clarify an instability of the ground state of the A chain against the lattice distortion that 
increases a strength (A) of a bond in each triangle. It relaxes the frustration and causes a remarkable 
gap enhancement; only a 6% increase of A causes the gap doubled from the fully- frustrated case 
(A = 1). The lowest excitation is revealed to be a kink-antikink bound state whose correlation 
length decreases drastically with A increase. The enhancement follows a power law, AiSgap ~ 
(A — 1) + 1.44(A — 1) 3 , which can be obtained from the exact result of the continuous model. This 
model describes a spin gap behavior of the delafossite YCUO2.5. 

65.40.Hq, 75.50.Ee, 75.60.Ch 



Cooperations between theory and experiment in the 
field of the low-dimensional quantum spin systems are of 
greater importance especially in these years, since vari- 
ous new compounds have been synthesized which realize 
theoretical models.EJ Much interest recently is focused es- 
pecially on the system with a spin gap. As tjapical exam- 
ples, we can consider the spin ladder modelj3 the bond- 
alternation model,El and the Majumdar-Ghosh model.cl 

The ground state of the A chain shares common prop- 
erties with models stated above. Its ground state is 
the pure singlet-dimer state with two-folded degenera- 
cies. The existence of the excitatioru energy gap above 
the ground state is rigorously proven.Q'Q However, strong 
frustration in this model plays an important role to cause 
unusual properties in the excited states, i.e., its excita- 
tion spectrum is dispersionless (zero-energy mode) .0 This 
mode contributes to the existence of the. additional low- 
temperature peak in thej-specific heat.ElB One of the au- 
thors (T.N.) and KuboJlS and Sen et alli^ clarified its 
elementary excitations and how they contribute to the 
thermodynamic properties. They found that the prop- 
erties of the excitation is governed by a kink and an 
antikink. The dispersionless aspect originates in a lo- 
calized kink, while an antikink moves as a free particle 
with an effective mass between two localized kinks. This 
picture is in a clear contrast with the excitation of the 
Majumdar-Ghosh model, whetc, both a kiek and an an- 
tikink are equivalently mobile.ta Sen et alOJ also pointed 
out the A chain can be a model for the recently discov- 
ered delafossite YCUO2.5, but the estimated spin gap is 
about half the one measured in the experiment. 

Strongly frustrated quantum spin systems have been 
attracting, since the interplay of the quantum effects and 
frustration may lead to an exotic ground state. Our ques- 
tion in this letter is that such a spin gap state caused by 
strong frustration is whether stable or not against a per- 
turbation. If relaxation of frustration occurs by a small 
perturbation, it may change the spin gap. A pairwise 
dimerization on the A bonds shown in Fig. |l| can be real- 
ized by the lattice distortion in this system, if the merit 
of the gap enhancement overcomes the distortion cost. 



Therefore we investigate the pairwise-dimerized A chain 
in this letter and clarify its gap behavior. The analyses 
are based upon numerical diagonalization, the variation 
and the continuous limit. First, we discuss the ground 
state and find a parameter region that can be realized in 
the material, and then we study the excitation. In due 
course, we will find a remarkable enhancement of the gap 
due to relaxation of frustration by a small increase of the 
pairwise bond strength (A); only a 6% increase of A causes 
the spin gap doubled, which is consistent with the gap 
measured in YCUO2.5. 



FIG. 1. Shape of the distorted A-chain. Double lines indi- 
cate the A bonds. 



We consider the following Hamiltonian as a model 
dimerized by lattice distortion. 



N 



(1) 



ri=l 



with 



AS^2n-l ■ S2n + 5^2n ' 5^2n+l + 5^2n-l ' S2n+1 (2) 



Here, N is the number of the triangles in the system, A 
is a parameter denoting the dimerization, and \S\ — 1/2. 
Figure ^ shows the depicted lattice. In this letter, we 
mainly consider the parameter A near the symmetric 
point A = 1. 

We numerically diagonalized the above Hamiltonian up 
to the systems with 28 spins {N = 14) under the periodic 
boundary conditions. The A dependence of the energy is 
shown in Fig. 0. 



1 



-9.5 I 

-10.0 





-10.5 




c 


-11.0 


u 






-11.5 


o 


-12.0 




-12.5 




-13.0 



k=0 R=+l 

k=0 R=-l 

k=7r R=-l 




0.8 



0.9 



1.0 



1.1 



1.2 



FIG. 2. The A dependence of the energy of the lowest state 
in each subspace denoted by k (wave number), and R (spin 
reversal symmetry) in the system with A'^ = 14 (28 spins). 
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FIG. 3. A dependence of the first excitation gap in A: = 
and k — TT sector. Sizes of the systems are N — S and — 14. 
For A < 1, gap in fc = tt becomes smaller than that in = 0. 



Tlie phase space is divided into a subspace character- 
ized by eigenvalues of the wave number (fc) and the spin 
reversal symmetry (i?) in order to reduce the memory 
use. From a technical reason, we only diagonalized the 
space with the wave number fc = and k — n, where the 
ground state and the relevant low-lying excited states are 
found when A = 1. It is noticed that the ground state at 
the fully-frustrated point, A = 1, is most unstable. 

The ground state in the case of A > 1 is trivial, i.e., the 
ground of the local Hamiltonian /i„ is the singlet dimer 
state on the A bond. Direct product of the dimers can 
span the whole system and be the exact ground state. 
Thus the total ground state energy is — 0.75AiV. Two 
excited states shown in the figure degenerate for A > 1, 
which was also observed for higher states. This degener- 
acy suggests that the excitation is local. 

On the other hand, the ground state for A < 1 is not 
trivially solved. It lies in the space of A; = 0, and the 
ground state energy gradually decrease with A go away 
from 1. A variational analysis and the second-order per- 
turbation analysis manifest that the energy difference be- 
gins with the order (1 — A)^. (Details are reported else- 
where.) Therefore the ground state is more stabilized by 
the perturbation with A > 1, which may be actually re- 
alized. In this letter, We only discuss the case A > 1 in 
detail, and leave the case A < 1 for another opportunity. 

Figure || shows the A dependence of the first excita- 
tion gap for both k = and fc = tt sector. Data of sizes 

= 8 and = 14 are compared to see the size depen- 
dence, as well. Degeneracy between "fc = 0" and "fc = tt" 
is solved for A < 1, and "fc = tt" becomes the first ex- 
citation. The gap enhancement is remarkably large for 
A > 1. Size dependence is also very weak for A away from 
1. At A = 1.06, the gap is about twice the one at A = 1. 



We consider the following open system to reveal the 
excited state. 

N 

^^opcn — hn + XS2N+I ■ S2N+2- (3) 

n=l 

We estimate the excitation energy by the following varia- 
tional analysis. The excited state we consider as the first 
approximation consists of (N — 1) singlet dimecs and two 
free spins. As was done in the previous paper ,E3 we name 
the free spins a kink and an antikink. The wave function 
of an antikink for an improved approximation may span 
beyond several spins since the free spin is not an eigen- 
state of the Hamiltonian. We have considered an antikink 
consists of 5 spins previously, but it merely increase the 
effective mass of an antikink and the shane,of the wave 
function are not significantly influenced. EJ'El 

We define a variational basis tpi so that an antikink is 
located at the ith triangle. 

= V-kink [4, 5] • • • [2i - 2, 2i - l]a2^ 

[2i+l,2i + 2]---[2N + 1,2N + 2], (4) 

where denotes a singlet dimer state connecting the 
ith and the jth site, namely [i, j] — aif3j — Piaj for ai{(3i) 
denoting up (down) spin located at the ith site. a2i is 
an antikink. The wave function of a kink located at the 
leftmost edge is known as, 

i/'kink = [ai(a2/33 - f32as.) + ^2(01/33 - /3ia3)]/V6. (5) 

This state is an eigenstatc of the local Hamiltonian hi, 
and therefore does not move. Its energy eigenvalue is 
A/4 — 1, therefore a kink contributes to the excitation 
by (A — 1). The singlet dimers that are not on the A- 
bonds, existing between a kink and an antikink, are not 
eigenstates of the local Hamiltonians, They also con- 
tribute to the excited energy. 

Variational bases are not orthogonal each other and 
satisfy the following relations. 
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(V'^IV'j) 



Si' 



ii;g + (A-l) + -(A-l)min(z,j) 

+ = Hi J 



(6) 
(7) 



Here, (5y is the Kronecker delta, and min(i, j) = i'lii < j. 
The ground state energy i?g = — 0.75A(-/V + 1). Since 
the matrix S given by eq. is the positive hermitian 
matrix, we can easily transform this variational problem 
with the trial function = ^ii^i into the eigen- 
value problem of the following Hamiltonian. 



(8) 



with 



4 



(9) 



where = min(«, j)^^. This can be solved numeri- 
cally for arbitrary N, say we show the result of = 200 
in this letter. Fig. ^ shows the wave function of the 
lowest eigenstate for A = 1.00, 1.001, 1.01. 
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FIG. 4. Variational wave function of the lowest eigenstate 
for A = 1.00, 1.001, and 1.01. Size of the system N = 200. n 
stands for the site number of the triangles. 



We find from this figure that an antikink drastically 
approaches a kink as A increase from 1. 

Fig. H shows the behavior of the gap enhancement 
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FIG. 5. Log-log plot of the gap enhancement for the vari- 
ational results of A'' = 200, and the numerical ones after ex- 
trapolated N ^ <x. We also plotted the exact result in the 
continuous limit of eq. ( |l2[ ) with the mass m — 1.21. 

We also plot the gap enhancement of the periodic sys- 
tem after the extrapolation of oo by using data 
with TV = 8, • • • , 14. Since the variational wave function 
of A = 1.01 only takes a finite value for n < 20, as shown 
in Fig. ^. the finite size treated in the numerical diagonal- 
ization is enough to extract the properties of infinite-size 
system for A > 1.01. Contrary to this, data for A < 1.01 
deviate from the variational ones, because the length of 
the bound state exceeds the finite size. The gap behavior 
qualitatively agrees with that of the periodic system, and 
is found to obey 

Sgap(A) - £;gap(l) ~ (A - 1) + 1.44 X (A - 1)2/3. (10) 

This behavior is what is expected by the following anal- 
ysis in the continuous limit of the model. 

We start by rewriting the variational bases and the rep- 
resentation of the Hamiltonian given by eq. into the 
momentum space. Since the wave length is defined to be 
sufficiently larger than the lattice constant in the limit, 
we consider the case of wave number k —f 0. We take into 
account up to the fc^ order for the diagonal part, and the 
leading order for the off-diagonal part. Then we obtain 
eigenvalue equations of the following Hamiltonian. 



ff(A).,. + (A-l)-J-|, + 5(A 



Da; 



(11) 



Here, a; is a distance between a kink and an antikink, m 
is an effective mass for an antikink, the constant term 



defined by £'gap(A) - £^gap(l) 



Eq{X) + £'gap(l) for Ea(X) is the ground state en- 
ergy. The third term of eq. ( [ll|) is a kinetic energy of 
an antikink, and the fourth term is a triangular potential 
expressing the energy loss by the singlet dimers between 
a kink and an antikink. Within the first approximation 
that an antikink is a free spin, the mass is 1 (m = 1). In 
the exact continuous model, on the other hand, structure 
of an antikink only renormalizes the mass as m(A), as far 
as the size of an antikink is negligible compared with the 
length of the bound state. 



3 



The eigenvalue equation H"^ = E"^ always gives bound 
states irrespective of A. This means that a kink and an 
antikink always form a bound state for A > 1. They col- 
lapse to a triplet state on a A bond in the limit A — + cxd. 
The wave function of the ground state of cq. (^l]) is given 
by the Airy function with the eigenvalue E: 

E-eo = iX-l) + l(^-^y\\-ir^'E (12) 

with E = 2.33816.0 We also plot this expression in 
Fig. ^ with the mass m — 1.21, the value estimated 
at A = 1,EJ though there is not a clear difference with 
the estimates using m = 1. This line coincide with the 
variational results for A close to 1, but it shifts upward 
as A increases. It suggests that the renormalized mass 
77i(A) > m(A = 1) = 1.21. An averaged length of the 
bound state is proportional to (A — 1)~^/^. 

We also calculated the magnetic susceptibility of fi- 
nite systems {N = 6 and 7) with periodic boundary con- 
ditions. All the eigenvalues were calculated by the nu- 
merical diagonalization. Therefore the numerically exact 
values were obtained. Figure ^ shows the temperature 
dependence of the susceptibility for A = 1.00, 1.06, 1.10. 

Peak position shifts toward the high temperature side, 
which is caused by the enhanced spin gap with increas- 
ing A. We consider the present data is enough to explain 
the thermodynamic limit for A = 1.06 and 1.10, since the 
length of the bound state is within the finite size. 
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FIG. 6. Uniform susceptibility per spin calculated for 
A = 1, 1.06 and 1.10 by the exact diagonalization of the finite 
systems with the periodic boundary conditions. The sizes of 
the systems are A'^ = 6(the dashed line) and 7(the continuous 
line). As A increase, size dependences become weak. 

In summary, we have investigated the A chain with 



the dimerized A-bond strength increased by a lattice dis- 
tortion. For the A > 1, the gap was found to obey 
(A - 1) + 1.44(A - 1)^/^; the first term is contributed 
by a kink and the latter is by an antikink. Only 6% 
change of the bond strength may cause the gap doubled. 
Therefore care must be paid for the comparison with the 
experimental results. Analyses on the case of A < 1 and 
the detailed derivation of the continuous limit will be re- 
ported elsewhere. 

Recently, Chitra et al investigated the 1-dimensional 
J1-J2 model with the bond alternation, and found that 
th&|gap behaves with 6^^^ with the alternation parameter 
6x3 We consider it is the same physics with the present 
model that a kink-antikink bound state is formed by the 
potential proportional to the 5 times distance between a 
kink and an antikink, which leads to the gap with expo- 
nent 2/3. 
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